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ABSTRACT 


An aircraft configuration for a deck launched intercept missicn was invesngated as part of an on- 
going wavender study by the Naval Postgraduate School and the NASA Ames Rescarch Center. The mission 
requirements for the carnicr-launched and recovered aircraft included Mach 6 cruisc out to a 1000 nautical mile 
combat radius and 20 minutes of combat followed by retum to the camer. A conical-flow wavender served as 
the starting point for the aircraft configuration. A hydrocarbon scramjct was integrated with the wavender 
body. The aft end of the waverider was faired to decrease the base area thereby reducing the transonic basc 
drag. A numencal optimization was then complcted to maximize the product of L/D (lift to drag ratio) and I,,, 
(specific impulsc). Variables for tie upumuzation included the cone shock angle (used to denve the conical flow 
waverider) and the gcometry of both the wavender body and the integrated propulsion system. The vehicle was 
constrained to a minimum volume of 3240 cubic fect, a maximum span of 60 feet and a fixed length of 60 foct. 
The integrated propulsion system was constrained to produce a minimum contraction rato of 12.0 and 
assurance that the cowl shock was within an acceptable distance of the inlct shoulder of the combustor. The 
optimum configuration met or exceeded all constraints. L/D comparisons were made between the integrated 
configuration (i.c., the subject of this study), pure Mach 6 optimized waveriders and histoncal trends. 
Additionally, model design, test media and test parameter sclection were studied for planned low speed wind 


and water tunnel tests as well as performance predictions for the planned wind tunnel tests. 
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I. INTRODUCTION 


The quest for speed has been a continuous theme in aeronautical history. The reasons 
for this quest have been many and varied; not the least of which has been the pure thrill 
and excitement of doing something in the present that could not be done in the past. 
Engineers, manufacturers, pilots and spectators have all shared in the excitement of 
increasing aircraft performance. Of course, the excitement of going faster does not, in 
itself, justify the cost. As noted by Kuchemann, there is a social aspect of aircraft design 
to consider [Ref. 1:p. 3]. Going faster must have a value to society which justifies the 
cost. Increases in speed have historically been worthwhile because they increased the 
convenience and decreased the time to travel a given distance, thereby increasing society's 
productivity. An increase in any figure-of-merit of aircraft performance must be balanced 
by other factors, such as productivity, cost, safety and environmental issues, in order to 
determine whether or not the increase is worthwhile. 

Kuchemann envisioned the peaceful future of our planet as a "global village" wherein 
which people could reach all of the world’s major population centers from any location in 
only a few hours. Kuchemann argued that this reduction in travel time was necessary to 
overcome the reluctance of people to travel great distances and that bringing about this 
"global village" by reducing travel times should be one of the social goals of aviation. 
Increasing the speed of aircraft may be the only way to achieve Kuchemann's vision. Frem 
our present status in aviation, hypersonic flight is the next logical step in the aircraft 
velocity spectrum. One possibility for efficient hypersonic flight is waverider 


configurations. [Ref. 1:pp. 448-451] 


To achieve his "global village", Kuchemann stated that transport aircraft of the future 
should be capable of traveling a “global range” (half the earth's circumference) in two 
hours. To this end, he defined a figure-of-merit by which to compare various aircraft 
types, as the fraction of the "global range" an aircraft was capable of traveling in two 
hours. Figure 1.1 shows a graphical representation of Kuchemann’s comparison R/Ry 
represents the fraction of the "global range" achievable in two hours by the different types 
of aircraft. Three general aircraft types are considered; swept, slender and waveriders 
The aircraft type is basically defined by its span to length ratio, s/, cruise Mach number, 
Mo, and body/wing shape, swept, slender or waverider. Three additional lines of interest 
subdividing the figure are the Mg=! line, the s//=0.2 line and the Bs//=I line Span to 
length of less than 0.2 leads to aircraft generally inadequate on the airfield (Ref 1 p. 449] 
Bs/=1 is where delta wings have a "nominally sonic leading edge." The figure shows that 
waveriders are best suited for the task of bringing the time for traveling “global range" 
down to a few hours. [Ref | pp. 448-451] 

Waverider technology has a wide range of possible applications, including, but not 
limited to, the well-publicized hypersonic transports of the future. The L/D performance 
of waveriders makes them a logical starting point for improving the performance of 
Maneuvering Reentry Vehicles (MaRV) and extending the range of Hypersonic Glide 
Vehicles (HGV) [Ref. 2:pp. 1-5]. Waverider characteristics also lend themselves to 
several space applications such as aero-assisted plane changes or orbit modifications and 
atmospheric long range cruising planetary probes [Ref. 3:p. 9]. Additionally, waveriders 
show promise as baseline configurations for single-stage-to-orbit vehicles. Their 
potentially superior L/D performance when operating at design Mach number compared to 
conventional wing/body configurations makes waveriders logical candidates for the next 


generation of hypersonic aircraft. 














Figure 1.1 Breguet Ranges of Various Aircraft Types. From Ref. 1. 


A. HISTORICAL PERSPECTIVE 

The history of waveriders apparently began in 1959 when Nonweiler introduced the 
idea of a three dimensional body derived from the flowfield behind a planar shock [Ref. 
4:pp. 521-528]. This first waverider had a delta planform and a caret shaped cross 
section. When flying at the design Mach number, the shock is attached to the leading 
edges preventing spanwise flow and spillage from the lower to upper surface. The high 
pressure flow behind the shock is captured between the shock and the lower surface of the 
body resulting in increased lift. Figure 1.2 shows Nonweiler's caret wing. 

In the 1960's, 1970's and 1980's, the design of waveriders was extended to include 
those derived from other known flowfields such as caused by right circular and elliptic 
cones. Figure 1.3 shows a conical flow waverider. In the late 1980's, Bowcutt and 


Anderson continued the progress in waverider development with viscous optimized 

















waveriders designed to maximize L/D. Figure 1.4 shows the Bowcutt/Anderson optimum 


result for a Mach 6 waverider design. [Ref. 6:pp. 15-19] 
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Figure 1.2 Nonweiler's Caret Wing. From Ref. 5. 


B. DESIGN MOTIVATION 
The configuration studied herein seeks to provide one possibility in the next logical 


step in waverider evolution; the design and optimization of a mission-specific waverider. 





As Schindel [Ref 7] points out, optimum waveriders are not feasible vehicles due to 
several practical considerations. These include volume requirements, the effects of base 
drag and the practicality of the proposed vehicle shape (construction and internal packing 
requirements). Additionally, an actual vehicle based on a waverider must be designed as 
an integrated engine/airframe unit. [Ref. 7:p. 1] 

In 1988 Vanhoy [Ref. 8] expanded the work of Bowcutt into the low-speed regime 
with wind tunnel tests of a Mach 6 viscous optimized waverider. Figure 1.5 shows the 
waverider used in Vanhoy's tests. In selecting the geometry for his tests, Vanhoy had to 
address such concerns as the ease of construction of the selected shape and whether the 
base thickness was adequate for use with a sting mount [Ref. 8:p. 8]. Such constraints 
represent only initial considerations in trying to create an actual vehicle from a waverider. 
The project described herein considers mission and performance constraints on the full 


scale vehicle as well as the requirements of a low-speed wind tunnel model. 
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Figure 1.3 Conical Flow Waverider. From Ref. 5. 
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II. MISSION DEFINITION AND ANALYSIS 


The carner-launched intercept mission was specifically chosen for the development of 
the waverider design discussed herein. Waveriders seem well suited to this type mission 
due to their inherent advantage in low drag at hypersonic speeds. Speed is an important 
factor in intercept missions because it translates into distance away from what the 
interceptor is protecting. Waveriders may also bring a better than average maneuverability 


at hypersonic speeds to the mission. 


A. MISSION SPECIFICS 

Building on the SR-71 and X-15 design, production and operations databases, a Mach 
number of six was chosen for the cruise speed as a logical, practical and technically 
feasible step up from current capabilities. Carrier basing imposes additional design 
constraints, including takeoff and landing speeds and weights, length and span and fuel 
type. It was decided that the more restrictive constraints of carrier operation during the 
initial design phase would lead to greater flexibility in relaxing constraints, if at a later date 
a similar land-based mission were desired. For example, design constraints such as length 
and weight isht be relaxed for a land-based waverider. Figure 2.1 depicts the proposed 
waverider intercept mission. Table 2.1 contains specific mission parameters chosen for 


this design. 








TABLE 2.1 
MISSION SUMMARY 


1000 nautical mile combat radius 
Cruise Mach 6 at Best Cruise Altitude (BCA 


~ = 


Twenty minutes on-station, airborne intercept, no dogfighting 













Loiter Twenty minutes at sea level at Best Loiter Speed (BLS) 
< 140 KTS 
Catapult launch and arrested recove 






— = 


Mach 6 Cruise Back Combat (28 min) 
Mach 6 Cruise Out 


Combat Radius 1006 miles ————_—__—_+, 





Figure 2.1 Mission Profile. 


B. AIRCRAFT SIZING 

An aircraft sizing study was completed for the selected mission parameters to be used 
herein for baseline comparison. The sizing was based on historical trends following 
Nicolai [Ref. 9] and Raymer [Ref. 10]. It is important to note that the sizing studies 
performed for this project are based on existing historical and current data bases. The lack 
of data on hypersonic aircraft makes predictions of aircraft size requirements and 


performance only somewhat better than an educated guess. It is, however, important in 





helping to validate conceptual results. Sizing began with weight, length and span 
limitations imposed by carrier basing. The sizing calculations are contained in Appendix 
A. A summary of preliminary sizing assumptions and results for the mission delineated 
above are given in Table 2.2. A comparison of preliminary sizing results with the 


optimized configuration is made in Section IV. 


TABLE 2.2 
PRELIMINARY SIZING SUMMARY 


St a ee OO NC ee et 









2670 square feet 
100,000 pounds 
Empty weight fraction 









C. ENGINE AND FUEL SELECTION 

For Mach 6 cruise, a supersonic combustion ramjet (scramjet) was selected as the 
appropriate propulsion system. Hypersonic aircraft require a more careful integration of 
the aircraft engine and body than conventional subsonic and supersonic aircraft due to 
their use of the aircraft body as compression and expansion surfaces for the engine. For 
the waveriders studied, the engine and body were designed from the beginning to work 
together as a single unit. Compression for the propulsion system is "mixed" in that it is 
accomplished by both internal and external means. External compression is accomplished 
by a series of oblique shocks which are created by ramps on the lower forebody. Internal 
compression is accomplished by a final oblique shock generated by the leading edge of the 
engine cowl. For flight below about Macl: 3 the scramjet is impractical and the use of 
more conventional engines becomes necessary. It is envisioned that two turbojets similar 


to J-58s could be used for takeoff, landing and low speed flight. For a 100,000 pound 
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aircraft this turbojet propulsion system would provide a thrust to weight ratio at takeoff of 


approximately 0.6. A longtime requirement of two engines for carrier-based aircraft has 
been standard practice in the U. S. Navy and is used in the mission analysis and sizing 
studies for this project. Table 2.3 contains a summary of engine sizing and performance. 
The two primary fuels considered for this design were hydrogen and ethylene. 
Ethylene, C2H4, was chosen as the fuel for this design for several reasons. One of the 
primary concerns for fuel selection was shipboard handling. Another major factor 
considered was the low density and therefore higher volume requirements of hydrogen. 
This creates problems for both storage space onboard ship and increased aircraft internal 
volume requirements to fly the same mission. Existing size constraints for carrier-based 
aircraft would limit the range of a hydrogen powered waverider. Finally, ethylene is very 
similar to JP-7; the fuel used by current J-58 turbojets. This allows the entire flight to be 


flown using a single fuel. 


TABLE 2.3 
ENGINE SIZING/PERFORMANCE SUMMARY 


Required T/W (takeo 
Required T/W (cruise 
Specific Impulse (hydrocarbon scramjet 
J58 Max takeoff thrust (sea level) 
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HI. WAVERIDER GENERATION AND OPTIMIZATION 


The waveriders analyzed for this project were generated using the Waverider Code 
developed by the Systems Analysis Branch at the NASA Ames Research Center. The 
Waverider Code and its subroutines are a subset of the Hypersonic Aircraft Vehicle 
Optimization Code (HAVOC), also created by the Systems Analysis Branch. The 
computer used was a Silicon Graphics Iris workstation. The design goal of this project 
was to achieve a hydrocarbon-scramjet-powered waverider optimized for mission 
performance, while keeping computer processing time to a minimum and checking the 
results enroute to solving for the optimum configuration. 

Due to the much increased computer processing time of the nose-to-tail flow solution 
by the Waverider Code when using hydrocarbon fuels (vice hydrogen), it was decided to 
begin by optimizing pure waveriders and then hydrogen-powered waveriders to generate a 
starting point for a hydrocarbon-powered waverider optimization. A non-optimized 
waverider had to be generated which met the optimization constraints to be used as a 
starting point for optimization due to difficulties experienced by the optimization routine 
when starting outside the defined constraints. This was accomplished by generating over 
one hundred waveriders to develop an understanding of how each variable affected the 
final waverider geometry and then finally moving toward a set of design parameters that 
produced a waverider which met the desired constraints. 

To save computer processing time and be able to more quickly determine an optimized 
waverider's suitability, it was decided to begin by optimizing a "pure" waverider (no 
engine, no base fairing) for maximum L/D. The result of this optimization would in tum 


be used as the initial geometry for a waverider with a base fairing, optimized for maximum 











L/D. The resulting optimized, faired waverider was used as the starting point for 
optimizing a hydrogen-scramjet-powered waverider. The objective function for this 
powered waverider was changed from L/D to the product of L/D and Isp. Using the 
product of L/D and Isp resulted in tradeoffs during the optimization process between 
aerodynamic and engine performance. Finally, the result of this optimization was used as 
the starting point for a hydrocarbon-scramjet-powered waverider optimized to maximize 


the product of L/D and Ixy. Figure 3.1 illustrates the design process used. 


13 


GENERATE NON-OPTIMIZED 
WAVERIDERS TO DEVELOP 
TRENDS 


DEFINE DESIRED 
PERFORMANCE, SIZE 
REQUIREMENTS AND 

CONSTRAINTS 


GENERATE NON-OPTIMIZED 
WAVERIDER WHICH MEETS 
CONSTRAINTS 


USE RESULTING NON- 
OPTIMIZED WAVERIDER AS 
START POINT FOR L/D 
OPTIMIZATION 


USE L/D OPTIMIZED 
WAVERIDER AS START 
POINT FOR L/D 
OPTIMIZATION OF FAIRED 
WAVERIDER 


INTEGRATE HYDROGEN 

POWERED SCRAMSJET ON 

RESULTING WAVERIDER 
AND USE AS START POINT 
FOR L/DxIsp OPTIMIZATION 


CHANGE FUEL ON 
RESULTING WAVERIDER TO 
ETHYLENE AND USE AS 
START POINT FOR L/DxIsp 
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Figure 3.1 Design Process. 








A. WAVERIDER GENERATION 

The Ames Waverider Code begins its computations by constructing a general cone- 
derived waverider from a user-defined flowfield as explained in Appendix B_ It departs 
from the resulting "pure" waverider with the addition of an integrated forebody ramp 
system, scramjet engine and nozzle. At the discretion of the designer, the lower surface is 
faired at the aft end to reduce the area of the base and thereby decrease the base drag at 
low supersonic and transonic speeds. The program then computes a nose to tail solution 
of the conical flowfield under the vehicle by solving the Taylor - Maccoll equation for 
conical flow. The code creates a data file which defines the geometry. This file 1s stored 
for later use to generate a three dimensional visualization of the body. 

The "nose to tail” solution of the conical flowfield includes a 2-D analysis of the inlet 
and nozzle flowfield and a quasi 1-D analysis inside the combustor. Both analyses 


consider real gas effects. Table 3.1 and 3.2 list highlights of each. 


TABLE 3.1 
COMBUSTOR FLOWFIELD ANALYSIS SUMMARY 


Augmented Pre-burning capabilit 

Heat balanced 

~ Ramp/Nozzle heat transfer rate input vs. Mach no. 

- Combustor skin friction coefficient input vs. Mach no. 





TABLE 3.2 
INLET/NOZZLE FLOWFIELD ANALYSIS SUMMARY 


Ramp/Cowl geometry input as a function of Mach no. 
Oblique shocks, isentropic waves and sliplines (contact surfaces) computed 
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B. OPTIMIZATION 

A numericai optimization was achieved using VMA Engineering's Design Optimization 
Tools (DOT). DOT is an optimization routine which varies user-selected parameters to 
find the local maximum or minimum value of a selected design objective function (i.e , L/D 
or L/D*I5p) while satisfying multiple user-defined constraints. It uses finite difference 
methods to calculate the gradients of the objective and constraint functions. DOT solves 
this nonlinear optimization problem iteratively. 

Specifically, DOT uses the Sequential Linear Programming Method to implement the 
Modified Method of Feasible Directions algorithm and find an optimum solution (in this 
case one which maximized L/D* Is, while still meeting all of the constraints). This method 
creates a Taylor series to approximate the objective and constraint functions and uses the 
resulting approximation for optimization, instead of the original nonlinear functions. This 
allows for easier calculation of the objective and constraint functions during the 
optimization. Additionally, the objective and constraint gradients can be taken directly 
from the linear Taylor Series expansion. DOT uses satisfaction of the Kuhn-Tucker 
Conditions to determine when it has reached a solution. Three additional considerations 
are checked to prevent unnecessary iterztions. The routine will end if a maximum number 
of iterations (defined by the user) is exceeded or if conflicting constraints prevent a 
feasible solution within 20 iterations. The third condition is a user defined tolerance which 
sets a minimum required change in the objective function after each iteration. If this 
condition is not met, asymptotic convergence is indicated and the optimization is 


terminated. [Ref. 1 l:pp. E-3 - E-36] 


C. OPTIMIZATION METHODOLOGY 
A specific, iterative methodology for the optimization process was selected to ensure 


consistency throughout the many optimization cycles. Figure 3.2 illustrates the selected 
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methodology. Each optimization began by defining the constraints, variables and objective 
function for the optimization and visually checking their assignment in the various code 
and data files. The affected files were then renamed according to run number and saved. 
This resulted in time savings and increased ease when trying to reproduce optimization 
results. The next step in the selected optimization process was to find a geometry which 
met the defined constraints, to act as an initial starting point for the optimization routine. 
Experience in how each variable effected the final geometry was used to generate these 
initial configurations. After running the optimization routine, the results were analyzed. 
This included checking to ensure that the waverider met all constraints, had an 
aesthetically pleasing geometry and delivered acceptable performance. If the optimization 
either failed or delivered unacceptable results, the variable or constraint responsible had to 


be determined and corrected and the process begun again. 
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Figure 3.2 Optimization Flowchart. 
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D. FIXED GEOMETRY PARAMETERS 

Several parameters were selected to remain fixed during the optimization process. 
These parameters were derived from mission and performance requirements and included 
the vehicle length, combustor design, fairing start point and equivalence ratio. 

1. Vehicle Length 

Vehicle length is constrained by aircraft carrier operations. A sixty foot length was 

selected as the maximum total vehicle length based on current carrier-based aircraft. 
Table 3.3 lists lengths of various carrier-based aircraft upon which the length selection was 


made [Ref. 12]. 


TABLE 3.3 
SIZE COMPARISON OF CARRIER-BASED AIRCRAFT 


AIRCRAFT MISSION LENGTH 
all weather attack | 53 ft O in 
62 ft 8 in 


E-2 Hawkeye early 80 ft 7 in 57 ft 7 in 
warning/control 


S-3 Viking 68 ft 8 in 53 ft 4 in 
F-4 Phantom II 38 ft lin 63 ft 0 in 
F-18 Hornet fighter/attack 37 ft 6in 56 ft 0 in 








2. Combustor 
The SCRAMJET combustor (powerplant) uses a three inch throat height and five 
foot combustor length. Combustor width was allowed to vary during the optimization, as 
discussed later. Maximum combustor efficiency was assumed to be 0.95 for this design. 
Combustor efficiency is the fraction of injected fuel to burned fuel. Tables 3.4 and 3.5 list 
locations of expansions and injectors inside the combustor. Figure 3.3 shows combustor 


geometry. Internal locations in the combustor are given by station numbers which are 


multiples of throat height (i-e., station number 1 is located 3 inches from the throat). 








TABLE 3.4 





















COMBUSTOR EXPANSION 
| Station | Distance from throat_| Area Ratio (Ai/A)) 
Er cee a | ee! 7, ee 2 ae 
eam ape (Pen oer F727 Cm 
Ag 20 in 

TABLE 3.5 
INJECTOR GEOMETRY 
be 25s deg a 





Figure 3.3 Fixed Combustor Geometry. 


20 





3. Equivalence Ratio 


The equivalence ratio is a measurement of how close to stoichiometric a reaction is, 


as defined by the following equation: 


= fuel | oxidizer 
© ( fuel | oxidizer) croicriomemne 


= 1 indicates a stoichiometric reaction, @ < | indicates a fuel lean (air rich) 
mixture and @ > | indicates a fuel rich (air lean) mixture. A stoichiometric reaction is one 
where all of the oxygen atoms in the oxidizer react chemically to appear in the products. 
[Ref. 13:pp. 63-69] A simplified example of the reaction seen with the combustion of 
ethylene is given by: 

C,H, + XO, + X(3.76)N, -> aCO, + bH,O +(3.76X)N, 
Reactants Products 
For a balanced equation, X, a and b must be defined such that: 


C,H, +30, + 3(3.76)N, > 2CO, + 2H,O +11.28N, 


The air to fuel (AFR) ratio necessary to provide stoichiometric combustion is 


determined by the following equation: 


AFR - New. - (343.7603) 


= 14.28moles air/ mole fuel 
Sud 1 


Ona weizt basis: 


AFR = 14.28(28.95) / 28 = 14. 76units air/ unit fuel 


An equivalence ratio of 0.6 was used for the optimization of the hydrocarbon- 
powered waverider. This value was chosen, after experimenting with different values, as 


the one which prevented thermal choking in the combustor and also resulted in the 


21 








optimizer driving the combustor width to an acceptable value. Thermal choking is defined 


as the supersonic flow in the combustor being reduced to Mach | by heat addition (Ref. 


14:pp. 77-85]. 
4. Fairing Start Point 
The parabolic fairing designed by Pessin (Ref. 15] was added to the lower aft end 


of the waverider (see Figure 3.4). The fairing was used to reduce the base area, thereby 





reducing the base drag for increased transonic and low supersonic performance. In doing 
this, however, several penalties were suffered. Among these were a reduction in internal 
volume and a loss of L/D performance. There were two basic considerations in selecting 
the fairing start point. A fairing start located farther forward from the base resulted in a 
sma'ter initial angle for the fairing, decreasing the effects of flow separation. The forward 
point was limited, however, by the requirement of the fairing to be concave upward. If the 
fairing began too far forward the result was a convex fairing, resulting in increased drag. 
Additionally, the cruise L/D performance went down as the length of the fairing was 
increased. This made a tradeoff between initial fairing angle and fairing length necessary. 
Since the optimization routine does not consider the increased problem caused by base 
area at transonic and low supersonic speeds, it was not possible to use the fairing start 
point as an optimization variable. Thus, experimentation with different values using the 
Waverider Code was required to determine an acceptable fairing start point. This 
experimentation led to a choice of start location at 85% body length or 51 feet from the 
aircraft nose. 
5. Assumptions 
Overall vehicle density was assumed to be 25 Ibs/ft3. A standard atmosphere was 


used for all calculations and the ratio of specific heats was assumed to be 1.4. 
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Figure 3.4 Pessin's Fairing Geometry. From Ref. 15. 


E. OPTIMIZATION CONSTRAINTS 

Constraints were placed on the optimization of the waverider in order to meet 
projected mission and engine requirements and to produce an aesthetically pleasing body 
shape. Specific constraints used during the optimization included: internal body volume, 
span to length ratio, contraction ratio, maximum leading edge temperature, level 
unaccelerated cruise, and a requirement that the cowl shock hit the combustor inlet 
shoulder. 

t. Volume 

Required volume was based on the preliminary mission analysis and sizing studies 

as discussed in Section II. The actual constraint used was a volume factor of 1.5% of the 
body length cubed. For a 60 foot long waverider this gives a minimum volume of 3240 
cubic feet. This results in an 81,000 pound aircraft at the assumed density of 25 Ibs/ft3. 

2. Span to Length Ratio 


A maximum S/L of 1.0 was chosen primarily to give a balanced appearance. 


Additionally, this leads to a reasonable maximum wingspan at the trailing edge of 60 feet 





for carrier based operations. Table 3.1 lists spans of various carrier-based aircraft for 
comparison. 
3. Contraction Ratio 
Contraction Ratio is a measure of exterr.al compression based on the size of the 
captured streamtube and throat height: 


CR= captured freestream area 
inlet area 


Contraction Ratio was constrained to a minimum of 12.0 for this optimization. 


Figure 3.5 shows the external compression shocks and contraction ratio geometry. 


Forebody 


Captured es 
Freestream > Te, 





Figure 3.5 Contraction Ratio Geometry. 


4. Unaccelerated Level Flight Requirements 
To achieve unaccelerated, level flight during the cruise phase of the mission, the 


forces acting on the aircraft must be balanced. Normally these forces are lift, drag, thrust 
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and weight. For unaccelerated flight, thrust (T) must equal drag (D). Similarly, for level 
flight, lift (L) must equal weight (W). At the speed and altitude considered for this 
mission, however, another force becomes significant. It is the centrifugal force generated 
by flying a flight path which accounts for the curvature of the earth. This centrifugal force 
acts in concert with the lift force to balance the aircraft's weight, thereby reducing the 
amount of lift required for level flight. Weight At Altitude (Wait) is defined as the aircraft 
weight minus the centrifugal force. Lift was required to equal Wa; during cruise. These 


requirements are summarized below and depicted in Figure 3.6. 


T=D 
Wy, =W-F 
L=W,, 


where: T=thrust 
D=drag 
W=weight 
F=centrifugal force 
L=lift 


Wait=weight at altitude 


5. Leading Edge Temperature 


Maximum leading edge temperature was constrained to 3250 °R. 





D T 
W 


Earth's curvature 
(exaggerated) 





Figure 3.6 Forces Acting During Level Unaccelerated Cruise. 


6. Shock on Shoulder 

During an early optimization run of a hydrogen powered waverider, the 
optimization process was requiring the cowl to be too far forward in order to increase the 
objective function. This resulted in increased compression, thereby increasing engine 
performance, by using reflected shocks inside the cowl prior to the throat. On an actual 
aircraft boundary layer effects behind these reflected shocks would probably lead to engine 
unstart. In order to prevent this process, a constraint was introduced into the optimization 
routine which required the cowl shock to fall on the inlet shoulder. This constraint was 
met by calculating a "Shock on Ramp" ratio. X, and Xg are the locations in percent body 
length of the cowl shock and ramp intercept and inlet shoulder, respectively. A tolerance 
of approximately 1 inch prior to the inlet shoulder was allowed in the optimization. For a 
body length of 60 feet, this leads to a Shock on Ramp ratio of 0.998. Figure 3.7 shows 


this "shock on ramp" geometry. 
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inlet shoulder 





Figure 3.7 Shock on Ramp Geometry. 


F. OPTIMIZATION VARIABLES 

Optimization variables included parameters which defined the waverider body and 
propulsion system geometry and the flowfield from which the waverider was derived. 

1. Generating Curve 

Input geometry variables include the first 6 coefficients which define the freestream 

trailing edge or generating curve (shown below). C7 and Cg were set to zero for this 
design. The generating curve is described by Figure 3.8. Table 3.6 lists coefficient limits 
for the optimization. 


X =C, tand, +C,Y + 10C,Y? + 100C,¥° + 1000C,Y‘ + 10,000C,Y* +C, cos(1.570796C,¥ / tand, ) 


where: 9, = shock angle of the conical flowfield 








Generating 
Curve 





Figure 3.8 Generating Curve Geometry (Aft View). 


TABLE 3.6 
GENERATING CURVE COEFFICIENTS 





2. Propulsion System Geometry 


A series of ramps on the underside of the aircraft is used to generate shocks which 
achieve the external compression of the air used by the scramjet. The geometry of these 
“compression ramps" is described by the angle and start locations of the ramps. Engine 


cowl geometry variables include the cowl start point and length. The initial nozzle angle 
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was also allowed to vary for the optimization. Table 3 7 lists the limits and initial values 
for these variables used during the optimization of the hydrocarbon configuration. Figure 


3.9 depicts propulsion geometry. 


TABLE 3.7 
PROPULSION SYSTEM GEOMETRY VARIABLES 


po tower limit | upper limit 
ramp istat | 0.27, | | 
framplangle [| 0 TT 
ramp 2 angle p00 8 
[shoulder location [ 0.602, | 0.69 | 
[cowlleadingedge [0.55 | 6s fT 
[nozzle start angle [200 | 40.0] 85 


angles in degrees 
locations in fraction of body length 











Body Length (60 f)_§ 


o—— lnitial Nozzle Angle 


| Lane Shoulder 
Cowl Leading Edge 


}-—— Cowl Length ————>| 





Figure 3.9 Propulsion System Geometry. 








3. Leading Edge Radius 
The leading edge radius of the waverider was allowed to vary between 0625 and 
1.0 inches. The initial value for optimization was set at 0.£ inches. A 0.5 inch leading 
edge radius was selected from previous waverider configurations as one which would 
result in a leading edge temperature of less than 3250 °R. 
4. Generating Flowfield Shock Angle 
The shock angle of the conical flowfield from which the waverider was derived was 
also allowed to vary during the optimization. The best L/D performance of several of 
Bowcutt's viscous optimized Mach 6 conical-flow derived waveriders occurred at a shock 
angle of 12 degrees (Ref. 16:p. 11]. From this information it was decided to let the shock 
angle vary near this value. From experience with prior waveriders it was noted that shock 
angles above 13 degrees led to poor L/D performance when the span to length was 
constrained to 1.0 and shock angles below approximately 11 degrees resulted in 
unacceptably low internal volume. Hence, the upper and lower limits were set at 13 and 
11 degrees, respectively. The initial value for the optimization calculations was 12.774 
degrees. 
5S. Fairing Start Angle 
The initial angle for the parabolic fairing was allowed to vary between 10 and 45 


degrees. The initial angle for optimization was set at 20 degrees. 


G. OBJECTIVE FUNCTION 

As discussed previously, the product of L/D and specific impulse was used as the 
objective function for optimization. The logic behind maximizing this value comes from 
the Bréguet Range Equation (for jets) shown below. Since the Bréguet Range accounts 


for propulsive efficiency, heat content of the fuel and aerodynamic performance (L/D), it 
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can be considered a “figure-of-merit for the whole aircraft" (Ref. |:pp. 6-7]. It can be 
seen from the equation that increases in either L/D or Isp lead to proportional increases in 


range for a given velocity and fuel. 


LW, 
Ren, Oy 
I tN, 


where: R = range 
H = calorific value of the fuel 
Np = propulsive efficiency 
WW, = segment weight fraction 


V = velocity 


Igp = Specific impulse 








IV. OPTIMIZATION RESULTS 


As discussed in Section III, the objective of this study was to optimize a hydrocarbon- 
powered waverider for the product of L/D and I,y, while meeting several size and 
performance constraints. A solution for an optimum hydrocarbon-powered waverider was 
found by the optimization routine which met all of the constraint functions. At the 
optimum solution, two of the seven constraint functions, Shock on Ramp and Contraction 
Ratio, were active. An active constraint is one in which the optimization routine is forced 
to "march" along the constraint function in order to increase the objective function and 
satisfy the constraint, rather than moving in the direction of the gradient. Fewer active 
constraints increases the likelihood of solving for the true maximum of the objective 
function in the local space being considered. The resulting flow properties, body 
geometry and propulsion system geometry of this optimum solution for a hydrocarbon- 


powered configuration are discussed below. 


A. FLOW PROPERTIES 
Table 4.1 contains free stream flow properties and flow properties at the cowl lip and 


combustor inlet for the optimum configuration’s calculated cruise altitude of 85,226 ft. 


TABLE 4.1 
OPTIMUM CONFIGURATION FLOWFIELD SUMMARY 


ee 
(psf) (R) ft/s) 








Free Stream | 6.00 | _0.3188_| 400.520 | 5886.2 


2.493 1474.563 | 4612.627 
Inlet 





















B. OPTIMUM BODY GEOMETRY 





Specific results from the optimized hydrocarbon-powered waverider body geometry 
are discussed individually below and are summarized in Table 4.3. As discussed in Section 
III, shock angle and generating curve coefficients were varied during tne optimization 
process. The values of these variables which led to the optimum configuration and the 
resulting body geometry are summarized below. Figure 4.1 depicts the resulting optimum 
configuration. 

1. Generating Shock Angle 

The angle of the conic shock which defines the flowfield from which the resulting 
optimum hydrocarbon-powered waverider is derived was selected as 12.9892 degrees by 
the optimization routine. 

2. Generating Curve 

The coefficients of the generating curve which defines the upper surface of the 
optimized hydrocarbon-powered waverider chosen by the optimization routine are listed in 
Table 4.2. It is noteworthy that all of the coefficients are positive. This leads to a simple 
curvature of the upper surface and would likely result in less expensive manufacturing 


costs for an actual aircraft. 


TABLE 4.2 
GENERATING CURVE COEFFICIENTS FOR OPTIMUM CONFIGURATION 


0.760425 | 0.049392 | 0.029621 | 0.047527 | 0.57519 | 0.014328 


3. Span 
The S/L of the optimized configuration is 0.9419. This was a direct result of the 
tradeoff between the contraction ratio constraint and trying to maximize L/D* Isp. Trends 


from the initial non-optimized waveriders indicated that, in general, L/D increased with 
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increased S/L. In order to achieve this increase in S/L, however, the waverider is 


generated lower (farther from the cone axis) in the conical flowfield. This leads to smaller 
contraction ratios and reduced engine performance. The resulting planform has a span of 
56.51 ft and a plan area of 2266.33 ft2. From this, Aspect Ratio is calculated below to be 
1.41 [Ref. 17:p. 158]. 





where: b = wing span 
S = plan area 
4. Internal Volume 
Internal volume is important because it determines the amount of fuel capable of 
being carried and, in turn, the aircraft's range. Volume Factor for the optimized 
configuration is 2.075 percent of the body length cubed. This leads to a volume of 
4481.04 ft3 and a weight of 112,025.93 lbs. 
5. Leading Edge Radius 
The leading edge radius of the optimized configuration is 0.4995 inches. This value 
is a result of the tradeoff during the optimization between leading edge temperature and 
drag. As the leading edge radius goes up, drag increases and temperature decreases. 
Similarly, as the leading edge gets sharper, the temperature increases and the drag 
decreases. It is desirable for a waverider to have as sharp an edge as possible in order to 


keep the shock attached to the leading edges. 
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Figure 4.1 Resulting Optimum Configuration. 
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TABLE 4.3 
OPTIMIZED CONFIGURATION BODY GEOMETRY 


Length Span Plan Area AR Volume Weight Lead. Edge 
Radius 


| oOf | 56.51 ft | 2266.33 Re | 1.41 | 4481.04 AF | 112,025.9 lbs | 0.4995 in 





C. OPTIMUM PROPULSION SYSTEM GEOMETRY 

The engine geometry parameters of the optimized hydrocarbon-powered configuration 
are summarized in Table 4.4. Combustor width of the optimum configuration is 10.238 ft. 
The forebody ramps, cow! and nozzle all share this width. The optimum configuration's 


propulsion system geometry is depicted in Figure 4.2. 


TABLE 4.4 
OPTIMIZED CONFIGURATION PROPULSION SYSTEM GEOMETRY 
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Figure 4.2 Optimized Configuration Propulsion System Geometry 
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D. OPTIMUM CONFIGURATION PERFORMANCE 

As discussed in the previous chapter, several constraints and an objective function 
were placed on the optimization process. All constraints were met by the resulting 
optimum hydrocarbon-powered configuration. Contraction Ratio and Shock on Shoulder 
Ratio were active at the optimum solution. Table 45 summarizes the optimized 


configuration's cruise performance as computed by the Waverider Code. 


TABLE 4.5 
OPTIMIZED CONFIGURATION PERFORMANCE SUMMARY 


CONSTRAINT LIMIT VALUE FOR OPTIMUM 
|___ContractionRatio [| 220 | 


Lift = Weight at Altitude (Wa) | Lift = 10654423 Ibs 
Way; = 106258.07 Ibs 














Drag = 27039.75 lbs 
Cruise D8 76 


E. COMPARISON WITH PRELIMINARY SIZING RESULTS 
Table 4.6 compares the optimized configuration with initial sizing results. The 


optimum configuration weight is for cruise. The preliminary sizing weights are at takeoff. 


TABLE 4.6 
COMPARISON WITH PRELIMINARY SIZING 








120,600 (first order 
| LENGTH TC OACSC‘C( ;CSdLSOOUOUOOOCCORTC 
|SPAN tt CT—“‘éORCUTS6 Gin 
| ASPECTRATIO. S| SECA 
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V. PERFORMANCE PREDICTION AND COMPARISONS 


Aerodynamic performance predictions were made using the VORLAX and Waverider 
codes at NASA Ames Research Center. The VORLAX code is a vortex lattice code that 
uses Polhamus'’ Leading Edge Suction Analogy. Briefly, Polhamus' theory accounts for 
flow reattachment inboard of the leading edge when calculating total lift [Ref. 8:p. 20]. 
This flow reattachment is caused by the normal force imparted by the vortex flow 
reenergizing the upper surface boundary layer. 

Lowspeed performance predictions (specifically C1.) are used in Section VI as a basis 
for wind tunnel model sizing and to predict forces measured in the tunnel. A comparison 
is made between predicted performance of the optimized hydrocarbon-powered 
configuration of this study, Vanhoy's lowspeed tested Mach 6 waverider (Ref. 8], one of 
Bowcutt and Anderson's Mach 6 viscous optimized waveriders [Ref. 16} and historical 


trends. 


A. COEFFICIENT OF LIFT 

The coefficient of lift, CL, is needed to predict forces seen in the wind tunnel. The 
first method used in predicting Cy, was low aspect ratio wing theory as given by the 
equation below. The resulting value multiplied by a (in radians) results in values for CL. 


versus a. 


AR = zit = 2.215 


Additionally, as noted by Pope [Ref. 18:p. 287], the Cy, curve for aspect ratios below 


three can be approximated by: 











dC, 
da 





= 0.008+0.018AR (per degree) 


Finally, a vortex lattice analysis was performed on the optimized hydrocarbon- 
powered configuration of this study by the Systems Analysis Branch, NASA Ames 
Research Center using the VORLAX code. One result of the VORLAX code is a 
prediction of Cy. at differing angles of attack and Mach numbers. Cy values were 
predicted by the VORLAX run for Mach numbers of 0.2 and 0.235 in anticipation of 
planned lowspeed wind tunnel tests. Angle of attack for these computations was varied 
from 0 to 20 degrees. 

The resulting Cy versus a values for both Low Aspect Ratio Theory and the 
VORLAX code are listed in Table 5.1. These results are also plotted in Figure 5.1. From 
the figure it can be seen that the slopes of the curves are almost identical. The values from 
the VORLAX code are slightly lower than those from Low Aspect Ratio Wing Theory. 
This is most likely a result of the slightly negative mean camber distribution of the 


waverider body as computed by the VORLAX code. 


TABLE 5.1 
VALUES FROM LOW AR THEORY AND VORLAX 


|__| Low ARTheory | Pope | 
LAOAP Mach 0.2 


-0.07546 
| 0.07145 | 


C 








0 Or 

0.07145 0.00421 
0.08414 
15.0 

20.0 


MIN 


~ 


| 10.0 | 0.3865 0.2618 0.24276 0.24295 
| 15.0 | 0.5798 0.3887 0.39534 0.39580 
| 20.0 | 0.7731 0.5156 0.53987 





0.53918 








| ——#®— LowAR Theory 
| —1+} — Mach 0.2 
——7 7 Mach 0.235 


——t— Pope 





3 
3 
€ 
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b= 
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Figure 5.1 Cy vs a. 


B. CRUISE PERFORMANCE 

Increased cruise performance makes waveriders generally better suited to hypersonic 
flight than other configurations. Of course, this advantage in cruise performance is based 
on the condition that the waverider be operating “on design"; the Mach number it is 
designed for and at zero angle of attack. When on design, waveiders have shown 
impressive results. Figure 5.2 shows a L/D comparison between three hypersonic 
configurations (Ref. 18:p. 28]. 

L/D at cruise for the optimum hydrocarbon powered configuration of this study is 
substantially lower than Bowcutt and Anderson's Mach 6 viscous optimized waveriders. 
As shown by Figure 5.3, the L/D of their best optimum Mach 6 waverider (0,=12°) is 
over 8.0 [Ref. 16:p. 28]. The lower L/D performance (4.17) of the waverider 
configuration of this study is directly attributed to the addition of the fairing and 


propulsion system. Table 5.2 summarizes this comparison. Cy for the configuration of 
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this study was computed at Mach six and zero angle of attack using the VORLAX code. 


The volumetric efficiency is given by: 


2 2 
“3 3 
qe I 
S, 2266.33 





where: V = internal volume 


Sp = plan area 
TABLE 5.2 
CRUISE PERFORMANCE COMPARISON 


[| _ Bowcut’s Optimum | Powered Optimum 
LD ee 











C. LOW SPEED PERFORMANCE COMPARISONS 

Comparisons of low speed results by Vanhoy [Ref. 8] indicated good agreement 
between Polhamus' Theory, a Mach 6 viscous optimized waverider, and a delta wing of 
similar length and span as shown in Figure 5.4 [Ref. 8:p. 41]. Similar performance should 
be expected by the configuration generated by this study. For comparison, the Cy, values 


computed by the VORLAX code have been overlaid on Figure 5.4. 
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Figure 5.2 Hypersonic Configuration Performance Comparison. From Ref. 19. 
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Figure 5.3 Results for a Series of Optimized Mach 6 Waveriders From the 
University of Maryland. From Ref. 16. 
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Figure 5.4 Cy, Comparisons. From Ref. 8. 








VI. TEST PREPARATION 


Low speed testing is important in determining the viability of any aircraft. It is critical 
for carrier based aircraft due to their low takeoff and approach speeds. Wind and water 


tunnel tests are planned for Summer 1993. 


A. WIND TUNNEL TESTS 

Low speed wind tunnel tests will be conducted using a 15 inch long aluminum model 
in the 32 inch by 45 inch wind tunnel at the Naval Postgraduate School. Low speed 
testing is intended to investigate the speed regime below approximately 150 knots; typical 
of carrier-based aircraft approach speeds. 

1. Wind Tunnel Description 

The wind tunnel planned for use is a closed circuit, single return, horizontal flow 

wind tunnel. Power is provided by a 100 horsepower electrical motor coupled to a three- 
bladed variable pitch fan via a four-speed truck transmission. The tunnel has a contraction 
ratio of 10:1 and a test section length of 48 inches. The maximum test section velocity is 
260 ft/sec with a nominal freestream turbulence level of 0.2%. Hinged glass windows on 
both sides and a window in the ceiling provide easy access and viewing. [Ref. 20:pp. 598- 


599} Figure 6.1 shows the planned installation of the waverider in the wind tunnel. 
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Figure 6.1 Planned Installation of Optimized Hydrocarbon-Powered Waverider in 
NPS 32x45 inch Wind Tunnel. After Ref. 20. 


2. Test Parameters 
Initial low speed testing will gather force and moment data from -10 to +20 degrees 
angle of attack. A sting mounted, six component internal strain gage balance will be used 
to measure the forces and moments experienced by the model in the wind tunnel. These 
measurements will yield lift, drag and side forces and pitching, rolling and yawing 


moments (Ref. 18:p. 126]. 


B. WATER TUNNEL TESTS 
Flow visualization tests will be performed in the Naval Postgraduate School Flow 


Visualization Water Tunnel Facility. This is a horizontal, continuous flow tunnel with a 15 
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inch wide, 20 inch high and 60 inch long test section. Glass construction allows easy 
viewing of the model during tests. Maximum water velocity in the test section is 
approximately one foot per second with a nominal RMS freestream turbulence level of less 
than one percent. Flow visualization is affected by dye injection at points of interest on 
the model. Six pressurized lines route the dyes from the dye reservoir to “he model. (Ref. 


21:pp 1-2] 


. PERFORATED INLET 

. DELIVERY PLENUM 

FLOW CONDITIONING ELEMENTS 
CONTRACTION SECTION 

OYE LINES 

. TEST SECTION 

MODEL SUPPORT 

. DISCHARGE PLENUM 
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Figure 6.2. The NPS Flow Visualization Water Tunnel Facility. From Ref. 21. 
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C. MODEL PREPARATION 
Model preparation included material selection, balance selection and model sizing for 
both the wind and water tunnel models. 
1. Wind Tunnel Model sizing 
Factors in choosing the model size included balance availability, wind tunnel test 
section size, and cost. The test section size of 32 by 45 inches limits the span of the 
model to approximately 18 inches. Appendix C contains model sizing calculations. 
2. Balance Selection 
A listing of available sting mounted, internal strain gage balances was obtained from 
the balance shop at NASA Ames Research Center. Several balance diameters were 
available. Due to the small model size required by the NPS wind tunnel, the 0.75 inch 
diameter balance was selected. An iterative process led to the selection of a 15 inch model 
length and a balance with a total normal force capacity of 50 pounds. At 20 degrees angle 
of attack, the model is predicted to generate a normal force of approximately 40 pounds 
(see Appendix C). 
3. Water Tunnel Model Sizing 
The major consideration in sizing the model for the flow visualization testing was 
flow blockage in the water tunnel test section. A limit of 5% flow blockage was used in 
sizing the water tunnel model. An 8 inch long model remains below this limit up to about 
22 degrees angle of attack (see Appendix C). 
4. Model Material and Construction 
Options considered for construction of the models included the computerized 
milling of aluminum, polyboard, or laminated mahogany models or the casting of epoxy 


resin models. Mahogany was ruled out for both wind and water tunnel tests due to its 


inability to hold the sharp leading edges required by the models. The radius of the leading 





edges on the wind tunnel model will be approximately 1 mil. Additionally, it would be 
difficult to seal the mahogany model for water tunnel testing. Polyboard was ruled out 
due to its lack of stiffness. Epoxy resin suits itself to making many copies of the same 
model since it is formed from a mold. The expense of using it to make just one model, 
however, was prohibitive. Aluminum was chosen for both models because of its rigidity, 
its ability to hold an edge and its relatively low cost. The specific alloy chosen was 7075 
aluminum, as recommended by the Ames machine shop. Table 6.1 contains material 


properties of 7075 aluminum alloy [Ref. 22:p. 245]. 


TABLE 6.1 
7075 ALUMINUM PROPERTIES 


Ultimate Tensile Young's 


Tensile eee Modulus 
Strength Streng 
at OD a 















Shear 
Modulus 
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VII. FOLLOW-ON RESEARCH 


The design of any new aircraft creates research opportunities in every engineering 
field. This is especially true if the aircraft being designed is on the frontier of technology. 
Hypersonic waveriders fit this description and work is required in every field to determine 
if they are viable candidates for the high speed aircraft of the future. Material and 
structural analysis, control system design, thermal control, and application studies are just 
a few of the areas that come to mind. For the waverider configuration of this study, the 
existence of a set of models makes intensive, hands on aerodynamic analysis possible in 


the near term. This should add dramatically to existing waverider data bases. 


A. PLANNED NEAR-TERM TESTING 

The next step in this joint, on-going study between the Naval Postgraduate School and 
the NASA Ames Research Center is the planned lowspeed wind tunnel and water tunnel 
tests for which the 15 inch and 8 inch models are currently being constructed. These tests 
are intended to gather force and moment data and flow visualization data, respectively. 
The results of both of these tests will undoubtedly lead to more opportunities for research. 

The initial water tunnel tests will likely show whether or not the position of the vortex 
flows around the leading edges of the waverider are stable. If they are not, testing to 
determine if they can be made stable can be studied. 

At 15 inches in length, the wind tunnel model is large enough to make the drilling of 


pressure taps at a later date relatively inexpensive. The material it is constructed of also 


makes testing in the higher subsonic and supersonic speed regimes possible. 











B. ADDITIONAL AREAS FOR RESEARCH 

Several specific questions arose during this study that apply directly to the waverider 
configuration designed herein. Each could be complete projects in themselves. Among 
these, modifications to the Waverider Code to append tails and/or a cockpit to the 
waverider is a possibility. Heat transfer, emissions, structural requirements, internal 


packing studies and refined mission performance predictions should be studied. 


Additionally, optimization of waverider configurations for other specific missions 


including single-stage-to-orbit, aero-assisted orbital transfer vehicles, missiles and various 
aircraft missions can be studied. Research in each of these areas, and many not listed, will 
serve to move waveriders closer to actual vehicles. Figure 7.1 depicts a conceptual 


aircraft based on the waverider configuration generated by this study. 











Figure 7.1. Hydrocarbon-Powered Waverider Aircraft Configuration. 
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Vill. CONCLUSION 


Arguably, the most valuable result of this research project has been the successful use 
of the Ames Waverider Code to generate a completely integrated engine/body waverider 
configuration optimized for a specific mission. The success of the Waverider Code in 
producing a viable powered waverider, with so many constraints placed on size, 
performance and appearance, attests to the code's robust nature and the programming 
talent of those who wrote it. 

It is concluded that an integrated engine/body aircraft configuration, capable of 
meeting specific mission requirements (a deck-launched interceptor in this case), can be 
derived from a conical-flow waverider. As expected, the aerodynamic performance of 
such a configuration, when compared with an optimum pure waverider, is substantially 
lowered by the addition of a propulsion system and fairing. Additional penalties will 
certainly be suffered with the addition of a cockpit, fin and control surfaces. 

This project has studied several considerations required to design, model and test an 
actual powered aircraft configuration based on a waverider. It is, however, only a first 
look at these requirements. It is recommended that studies in the optimization of 
waveriders for application to other real-world missions be undertaken. These studies 
should be of integrated engine/body designs with a view towards optimizing a specific 
desired mission performance. Only this type of research will answer whether or not 


waverider configuratiuns can be viable vehicles. 
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APPENDIX A - PRELIMINARY AIRCRAFT SIZING 


A preliminary sizing study was done following Nicolai (Ref. 9] and Raymer (Ref. 10]. 
Preliminary sizing was based on historical trends and current data bases and was used for 
comparison with the optimized hydrocarbon configuration of this study. The basic result 
of preliminary sizing is to estimate gross takeoff weight and fuel weight required for the 


mission. 


A. WING GEOMETRY 
Wing geometry is predicted for the given mission based on historical trends. 
1. Aspect Ratio 
One would expect that a delta-type planform waverider would have a low aspect 
ratio. Based on historical trends for jet fighters (non-dogfighters), the aspect ratio is 
estimated at 1.35 [Ref. 10:p. 53]. 
2. Wing Area 
Based on the predicted aspect ratio and assuming the span is equal to 60 feet 


(S/Lmax=!.0), the wing area, S, is equal to 2670 square feet. 


2 2 
Sa 67072 
AR 1.35 


B. THRUST-TO-WEIGHT RATIO 
1. Takeoff 


A statistical estimation of T/W for a jet fighter (non-dogfighter) gives a typical 


required installed T/W of 0.6 [Ref. 10:p.79]. Using two J-58 turbojets, each with 30,000 














pounds of thrust at sea level [Ref. 9:p. 14-6], fills this requirement for a 100,000 pound 
aircraft. 
2. Cruise 
During level, unaccelerated cruise, the T/W must equal the inverse of L/D (Ref. 
10:p. 81]. Cruise L/D was assumed to be 4.0 for this calculation. 


(7/W] ce = soe, Ss 


[L/D] cu 


C. TAKEOFF WEIGHT ESTIMATION 

Takeoff weight is estimated using an iterative process based on an initial guess of 
takeoff weight, crew and payload weight, and calculated fuel and empty weight fractions 
as given below [Ref. 10:p. 12]. Fuel fraction assumes a 6% allowance for reserve and 
trapped fuel. (Ref. 10:p. 23} 


W. = Woew + W sayioad 
° 1-(W,/W,)-(W./,) 





i = 1.06(1- 7 ) 
where: Wo = takeoff weight 
Werew = Crew weight 
Wpayload = payload weight 
We = fuel weight 
We = empty weight 


W, = mission segment fuel fractions as calculated below 








1. Initial Takeoff Weight Estimate 
The initial weight guess for this aircraft was chosen to be 100,000 pounds based on 


the maximum takeoff gross weight for existing catapults from Figure A.1 [Ref. 10:p. 90} 


CATAPULT 
END SPEED CATAPULT TYPE 
TS) Cll C7 C13 
150 





40 60 
MAXIMUM TOGW, 10 Ibs 








Figure A.1 Catapult End Speeds. From Ref. 10. 


2. Empty Weight Fraction 
Empty-weight fraction is estimated to be 0.5 for initial sizing. This is based on 
historical trends and also on catapult takeoff and arrested recovery weight limitations. 
Maximum catapult takeoff weight is approximately 100,000 lbs. From Figure A.2, 
maximum arrested recovery weight is approximately 50,000 Ibs [Ref. 10:p. 92]. These 


maximums are used as a conservative estimate for initial sizing. The use of advanced 








composites in construction of the aircraft could result in a 16% lighter empty weight (Ref. 
9:p. 5-2]. This results in an empty weight fraction of. 


W,_ (50,000)(1-0.16) _ 5 4, 
W, 100,000 


MARK 7 ARRESTING GEAR 


MOD 3 


MOD2 
MOD 1 


20 40 60 80 


MAXIMUM LANDING WEIGHT, 10° Ibs 





Figure A.2 Arresting Gear Weight Limits. From Ref. 10. 


3. Mission Segment Weight Fractions 
a. Warm-up and Takeoff 
Weight fraction for engine start, taxi, and takeoff is estimated from historical 
trends [Ref. 10:p. 104]. 


W /Wo = 0.97 
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b. Climb and Accelerate 


For supersonic aircraft the climb weight fraction with beginning Mach of 0.1 Is 


approximated by: 


a = 0.991-0 007M -0.01M* 


i-| 


where: M = ending Mach number 
The beginning Mach number for this design will be based on the catapult end 
speed for the aircraft weight. From Figure A.1, this is approximately 130 kts or M = .197. 
For a cruise Mach number of 6.0 the weight fraction for acceleration from M = 0.1 must 
be divided by the weight fraction for acceleration from M = 0.1 to M = 0.197. [Ref. 10:p. 
105] 


W, __ 0.991—0.007(6.0)—0.01(36) _ 


W, 


= 0.595 
0.991 - 0.007(.197) — 0.01(.039) 


c. Cruise 


The Brequet range equation is used for cruise: 
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This leads to cruise weight fractions of: 


-RC 
FCLID) 


W, 
—_=e 
W, 
where: R = range 
C = specific fuel consumption 
V = velocity 


L/D = lift-to-drag ratio 


(1) Range. Range for each cruise segment (out and back) is 1000 nm. 

(2) Specific Fuel Consumption. From Figure A.3, an initial estimate of Isp 
equal to 1300 seconds is made for a hydrocarbon-powered scramjet. Specific fuel 
consumption, the inverse of specific impulse, is equal to 7.7x10-9/sec. 

(3) Velocity. Velocity is calculated assuming a standard atmosphere at 80,000 


feet. 


v= Ma 


where: M = cruise Mach number 
a = speed of sound 
At 80,000 feet the speed of sound is 968.08 ft/sec. For a Mach number of 6.0 
this gives a velocity of 5808.48 ft/sec. 
(4) L/D. L/D is assumed to be 4.0 for initial sizing calculations. This is a 


conservative assumption based on trends developed during early waverider configurations. 
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Figure A.3 Specific Impulse Trends. 


(5) Cruise Out 


- (6016000 ft)( 7.7*10 “ see) 
(5808.48 ft/sec)(4.0) = 0.818 


(6) Cruise Back 


cd. Combat 
Combat fuel requirements are a function of engine performance and combat 


duration [Ref. 10:p. 106]. For the specified mission, the combat duration is 20 minutes. 





It is assumed for these calculations that combat will be flown at low Mach numbers using 


the turbojets for propulsion. T/W for this segment is estimated to be 0.5. 


7 =1-C(T/W)(d) = 1- 1.72 x 10 *(.5)(1200) = 0.897 
where: C = specific fuel consumption (1/sec) 
T/W = combat thrust to weight ratio 
d = combat duration (sec) 
e. Loiter 
Typical requirements for aircraft include the capability of 20 minutes of loiter 
before landing. This allows for such eventualities as missed approaches and fouled 
runways/decks. Specific fuel consumption for this mission segment is based on the J-58 
turbojets and is equal to 1.72x10-4/sec. Loiter for jet-powered aircraft is flown at 


maximum L/D. Maximum L/D is estimated to be 8.0 for these calculations. 


Ww. ~EC ~(1200sec)(1.7210*) 
wot =e ae = 0.975 
§ 
where: E = loiter time 


Sf. Descent for Landing 


Descent for landing is estimated from historical trends (Ref. 10:p. 107]. 
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g. Landing and Taxi Back 
Landing is estimated from historical trends [Ref. 10:p. 107}. 


Mh 992 
W. 


3. Fuel Fraction 
The fuel fraction is calculated from the total mission weight fraction given by the 


product of the mission segment weight fractions calculated above. 


8 — (.97)(.595)(.818)(.897)(.818)(.975)(.990)(.992) = 0.332 


= 1.06(1 — 0.332) = 0.708 


R\S 31 


4. Resulting First Order Estimation of Takeoff Weight 


W, = W, + W, load Wa tW ow, 


pay: 


= 200 + 400 + (.708)(100,000)+.42( 100,000) = 120, 600/bs 
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APPENDIX B - WAVERIDER CONSTRUCTION 


Waveriders are generated in a reverse fashion from most aerodynamic bodies in that 
the flowfield from which the waverider is derived is known. Waveriders can be generated 
from the flowfields caused by shapes other than themselves (wedges, cones, etc.) The 
waveriders of this study were generated from the flowfields caused by axisymmetric 
circular cones at zero angle of attack. The explanation of waverider construction 
contained herein will use this type of waverider as an example. The process used for the 


waveniders of this project is straightforward and fairly standard. 


STEP 1. SOLVE THE FLOWFIELD 
The first step in creating a conical flow waverider is to solve for the flowfield 
generated by a right circular cone of known semi-vertex angle. The Taylor Maccoll 


equation (Ref. 16:p. 7] is used to generate the solution for this flowfield: 


(y-1)/2[t-V? - (av, 1d0)' |[2V, + (dV, /dO)cot 0+ d?V, 1d@ - (dV, 1d) + (dV, /dO)(d°V, 1d0*)] = 0 


where: V, = nondimensional component of flow velocity along a conical ray 
© = ray angle from axis 


Y = ratio of specific heats 


STEP 2. GENERATE FREESTREAM UPPER SURFACE 
The generating curve (see Section III) is drawn longitudinally through the cone 
flowfield to "carve" out the waverider's freestream upper surface. Figure B.1 shows a 


general cone-derived waverider, the generating cone and the conic shock. 








CONE VAHCH GENERATES 






freestream surface 
compression surtace 


Figure B.1 Construction of General Cone-Derived Waverider. From Ref. 16. 


STEP 3. LOWER COMPRESSION SURFACE 
Conic streamlines are traced back from the intersection of the shock and the upper 


surface to generate the lower, compression surface of the waverider. 
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APPENDIX C - MODEL SIZING CALCULATIONS 


A. WIND TUNNEL MODEL 

The two major factors considered in sizing the wind tunnel model were the test section 
size and normal force measuring capability of the sting mounted balance. The test section 
size of the NPS tunnel (32 x 45 incies) limits the span of models to approximately 18 
inches. The small size of such a model indicated the need to use a 0.75 inch diameter 
balance (the smallest available from Ames). The 0.75 inch diameter balances were 
available with normal force capabilities of 50 and 100 pounds. A general rule of thumb for 
force measurements is to use approximately 75 percent of the balance capacity at the 
highest force measured. This will generally result in the best resolution of measured 
forces. The goal of this sizing analysis was to find the largest model possible which would 
generate a normal force of 75 percent while at the maximum angle of attack planned for 
the tests. 

The normal force expected in the tunnel is a function of the dynamic pressure, q, the 
coefficient of lift, Cy, and the size of the model. For best resolution, this force should be 
37.5 or 75 pounds for the 50 and 100 pound balances, respectively. The maximum angle 
of attack for the initial tests will be 20 degrees. The predicted Cy value at 20 degrees 
angle of attack was calculated by the VORLAX code at NASA Ames as 0.54. 

For a standard atmosphere at sea level, the speed of sound is 1116.4 ft/sec and the 
density is 2.3769x1073 slugs/ft3. The planned test velocity is 140 kts or 236.3 ft/sec. 


From this, the dynamic pressure is calculated as follows: 





q= lad = (.5)(2.3769 x 10° )(236.3)? = 66.36Ibs / ft’ 


where: p = density (slugs/ft3) 
v = velocity (ft/sec) 


From the optimized hydrocarbon-powered configuration, the plan area, s, is given by: 


s = (.674)(1)? (.9416) 


The lift force can be calculated by: 


L=C,qs 


where: Cy = coefficient of lift 
q = dynamic pressure 
s = plan area 
The 100 pound balance was ruled out due to the inability of an 18 inch model (the 
largest possible in the tunnel) to generate 75 pounds of normal force (lift) at 20 degrees 


angle of attack as shown below: 


L = (.54)(66.36)(1.428) = 51.2 Ibs 


The model length necessary to generate 37.5 pounds of force can be found in a reverse 


I= ee = 1.28 ft = 15.4 in 
(.54)(66.36)(.674)(. 9416) 


fashion: 








B. WATER TUNNEL MODEL 

The major factor in sizing the water tunnel model is flow blockage by the model in the 
water tunnel test section. A general rule of thumb is to limit this flow blockage to five 
percent of the test section area. The test section area of the NPS water tunnel facility is 
300 in2. This gives a desired maximum blockage of 15 in2. The blockage is a function of 
the model length and angle of attack. For this series of tests the maximum angle of attack 
of interest is approximately 20 degrees. The maximum length of the model can be 
calculated by knowing the plan area which generates the maximum blockage at 20 degrees 
angle of attack and is calculated below. 


The maximum blockage is given by: 





(.05)(300in’) = 15in? 
From this, the plan area of the model follows as: 


15in” 


sin(20°) = 43.86in’ 
Sin 


The relations between length, span and plan area are known from the optimized 








hydrocarbon-powered configuration and are: 


plan area = (.674)(length)(span) 


span = (.9416)(length) 


From this, the model length is computed as: 


length = Regie eee = 8.3in 
(.674)(.9416) 
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From the results above, a length of eight inches was selected. This length will result in 


a maximum angle of attack of approximately 22 degrees at five percent flow blockage. 
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